Energetics of PCMDI/CMIP3 Climate 
Models: Energy Budget and Meridional 
Enthalpy Transport 

V. Lucarini, Dept. of Meteorology & Dept. of Mathematics 
University of Reading, Reading, RG6 6BB, UK 
Email: v.lucarini@reading.ac.uk 
F. Ragone, Meteorologisches Institut, KlimaCampus 
University of Hamburg, Hamburg, Germany 

October 15, 2011 

Abstract 

We analyze the PCMDI/CMIP3 simulations performed by climate 
models (CMs) using pre-industrial and SRESA1B scenarios. Rela- 
tively large biases are present for most CMs when global energy bud- 
gets and when the atmospheric, oceanic, and land budgets are consid- 
ered. Apparently, the biases do not result from transient effects, but 
depend on the imperfect closure of the energy cycle in the fluid com- 
ponents and on inconsistencies over land. Therefore, the planetary 
emission temperature is underestimated. This may explain the CMs' 
cold bias. In the pre-industrial scenario, CMs agree on the location in 
the mid-latitudes of the peaks of the meridional atmospheric enthalpy 
transport, while large discrepancies exist on the intensity. Disagree- 
ments on the location and intensity of the oceanic transport peaks are 
serious. With increased CO2 concentration, a small poleward shift of 
the peak and an increase in the intensity of the atmospheric transport 
of up to 10% are detected in both hemispheres. Instead, most CMs 
feature a decrease in the oceanic transport intensity in the northern 
hemisphere and an equatorward shift of the peak in both hemispheres. 
The Bjerkens compensation mechanism is active both on climatolog- 
ical and interannual time scales. The peak of the total meridional 
transport is typically around 35° in both hemispheres and scenarios, 
whereas disagreements on the intensity are relevant. With increased 
CO2 concentration, the total transport increases by up to 10%, thus 
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contributing to polar amplification. Advances in the representation of 
physical processes are definitely needed for providing a self-consistent 
representation of climate as a non-equilibrium thermodynamical sys- 
tem. 

1 Introduction 

The investigation of the global structural properties plays a central role for 
the provision of a unifying picture of the climate system (CS). Such an en- 
deavor is of fundamental importance for understanding climate variability 
and climate change on a large variety of scales, which encompass major pale- 
oclimatic shifts, almost regularly repeated events such as ice ages, as well as 
the ongoing and future anthropogenic climate change, as envisioned by the 
scientific programme proposed in the landmark book by Saltzmann [47]. 

This effort has significant relevance also in the context of the ever- increasing 
attention paid by the scientific community to the quest for reliable metrics 
to be used for the validation of climate models (CMs) of various degrees 
of complexity, as explicitly requested by the 4th Assessment Report of the 
Intergovernmental Panel on Climate Change (IPCC4AR) [51] and for the def- 
inition of strategies aimed at the radical improvement of their performances, 
beyond incremental advances often obtained at the price of large increases in 
requested computing power. See the discussions in, e.g., [27, 16, 35, 36, 28]. 
EU projects such as PRUDENCE, ENSEMBLES, DEMETER, US initiatives 
such as the Phase 3 of the Third Coupled Model Intercomparison Project 
developed within the Program for Climate Model Diagnosis and Intercom- 
parison (PCMDI/CMIP3) and innovative volunteering activities as those or- 
ganized by the UK climateprediction.net project have fueled great efforts 
aimed at auditing CMs' performances. The pursuit of a "quantum leap" in 
climate modeling - which definitely requires new scientific ideas rather than 
just faster supercomputers - is becoming a key issue in the climate community 
[50]. 

In a modern perspective, the climate can be seen as a complex, non- 
equilibrium system, which generates entropy by irreversible processes, trans- 
forms moist static energy into mechanical energy [25, 26] as if it were a heat 
engine [19], and, when the external and internal parameters have fixed val- 
ues, achieves a steady state by balancing the input and output of energy 
and entropy with the surrounding environment [54, 44, 43]. The tools of 
phenomenological non-equilibrium thermodynamics [46, 6] seem very well 
suited in defining a new point of view for the analysis of the CS [21, 29], 
for understanding its variability and its large-scale processes, including the 
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atmosphere-ocean coupling and the hydrological cycle. Moreover, such ap- 
proach seems valuable for understanding the mechanisms involved in climate 
phase transitions observed at the so-called tipping points, i.e., conditions 
under which catastrophes may occur for small variations in the boundary 
conditions or in the internal parameters of the system [24]. Recent results 
suggest that great progresses in this direction are definitely within reach 
[37, 38]. 

The energy driving the CS comes from the Sun in form of shortwave 
radiation. Part of the radiation is scattered elastically back to space, with 
an efficiency which depends critically on the surface irradiated the incoming 
light. Clouds, ice- and snow-covered and barren surface have in general the 
ability to scatter back to space relevant portions of the incoming radiation, 
while the process is much less efficient for, e.g., forested lands and especially 
the ocean surface. Albedo is commonly referred to as the fraction of the 
scattered with respect to the incoming radiation, and is in general a two- 
dimensional time-dependent variable. The rest of the incoming radiation is 
absorbed by the fluid and solid components of the CS. The absorption process 
is not homogeneus in space and time, because of variability in the amount 
of incoming energy and because the interaction of radiation with matter de- 
pends critically on the nature and density of the irradiated medium. Most of 
the net incoming solar radiation is absorbed at surface and is redistributed 
throughtout the atmospheric column by a variety of means such as con- 
vection, latent and sensible heat fluxes, and longwave radiative exchanges 
within the atmosphere and between the atmosphere and the underlying sur- 
face. Organized and turbulent atmospheric motions and the hydrological 
cycle (including the formation and dissipation of clouds) are the direct result 
of this process of energy redistribution. Oceanic motions are also generated, 
basically as an outcome of the variable exchange of heat, water and mo- 
mentum between the water masses and the overlying atmosphere. Energy 
is irradiated back to space in the form of longwave radiation approximately 
thermalized at temperatures typical of the CS, so that on the global aver- 
age the net incoming solar radiation is approximately compensated by the 
outgoing longwave radiation. Since clouds have a great impact in both de- 
termining the amount of scattered solar radiation - because of their typically 
high albedo - and of outgoing longwave radiation - because of their typically 
high emissivity - they have a crucial role in modulating the local and global 
energy budget. On a planetary scale, large scale atmospheric and oceanic 
motions allow for the reduction of the temperature difference between high 
and low latitudes in both hemispheres with respect to what would be en- 
forced by purely local radiative and convective processes. The low and high 
latitude regions differ largely in terms of net incoming radiation, because of 
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astronomical factors and local albedo properties, and feature opposite sign 
of the net budget of the radiative flux at the top of the atmosphere (TOA). 
In general, we understand that the so-called differential heating, both in the 
small and in the large scales, is the basic ingredient of the non-equilibrium 
nature of the CS [26, 44]. The described processes have relevant impacts 
also in terms of the second law of thermodynamics, as in the CS entropy 
production is inextricably related to the irreversible heat transfer from high 
to low temperature regions [14, 43, 37, 38]. 

Therefore, it seems of the greatest importance to devote attention to the 
analysis of the consistency of the state-of-the-art CMs in the representation 
of the energy and entropy budget and transport in the CS as a whole as well 
as in its main subdomains. While the analysis of entropy budget will be the 
subject of a separate publication, we hereby wish to focus on the analysis 
of energetics of the CS in selected simulations of the CMs included in the 
IPCC4AR, by taking advantage of the immense data repository provided by 
the PCMDI/CMIP3 project. 

In this work we focus on one side on the behavior of the CMs under steady 
state conditions, and thus take advantage of the data produced with the pre- 
industrial scenario, where CMs are integrated for several hundreds of years 
with fixed atmospheric composition with 280 ppm CO2 concentration. On 
the other side we aim at understanding how different is the CMs' response 
to the CO2 radiative forcing impacts. For this latter purpose, we have con- 
sidered the CMs' outputs produced under the Special Report on Emission 
Scenario A1B (SRESA1B) scenario, which foresees a stable CO2 concentra- 
tion of 720 ppm after 2100 following a prolonged increase after 2000. The 
initial conditions of the SRESA1B runs are provided by the final state of 
the XX century runs, which are in turn initialized in the second half of the 
XIX century using the pre-industrial runs state as initial conditions. Since 
we wish to avoid the period of direct forcing and focus on the time frame 
when CMs start to relax towards the steady state, we analyze data referring 
to the XXII and XXIII century. 

1.1 Global Budgets in a Steady State 

Following and expanding previous analyses [66, 65, 63, 59, 61], we first per- 
form a thorough audit of the CMs in their representation of the energy budget 
of the whole CS, and, separately, of the atmosphere, of the ocean, and of the 
land, in order to test how consistently steady state conditions are realized in 
the case of pre-industrial simulations. Moreover, by comparing the obtained 
results with the energy budgets obtained from the runs performed under 
SRESA1B conditions, we understand how the CMs describe the CS response 
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to time- dependent CO2 radiative forcing. 

As discussed above, basic physical principles suggest that under generic 
boundary conditions the ultimately realised steady state obeys general zero- 
sum properties for energy fluxes [46, 44]. Therefore, when statistical sta- 
tionarity is eventually obtained after transients have died out, a physically 
consistent - even if not necessarily realistic in the representation of each and 
every climatic feature - CM should feature a vanishing (when long time aver- 
ages are considered) energy budget as a whole and for any of its subdomains. 
If this is not the case, the energy bias of the CM has to be attributed to 
unphysical energy sources or sinks. We remark that this property has noth- 
ing to do with the so-called operation of model tuning: a CM with a bad 
parametric tuning probably features unaccurate or even unrealistic climate 
properties, but, if the CM is energetically consistent and steady state is at- 
tained, the energy budget for any subdomain should be vanishing on the 
average. 

The first panel of Fig. (1) contrasts the idealized behavior of a CM whose 
energy exchange processes are consistent (solid line) with that of a biased CM 
(dashed line). In the former case, in the initial steady state conditions (black 
line) the TOA average energy budget vanishes, whereas in the latter a spu- 
rious bias (Bias 1 in the figure) is present. When the parameters of the 
system are changed with time so that forcing is applied (red line), both mod- 
els foresee a non vanishing TOA budget (positive in this representation, as 
in global warming conditions). When the parameters are held fixed to newly 
established value, the system relaxes toward the new steady state within a a 
time comparable to a few units of the slowest of the internal time scales (the 
oceanic one, in usual climate conditions). But whereas the energy-consistent 
model readjusts towards a zero TOA energy balance, the other model's TOA 
energy budget will readjust towards a non- vanishing bias, generically differ- 
ent (and so climate-dependent) from the initial one. Note that, whereas in 
the case of a good model the net total heating resulting from the forcing 
corresponds to integral of the energy budget, in the case of imperfect model 
with potentially climate-dependent bias the net heating of the system is in 
principle ill-defined. Conceptually analogous cartoons actually apply not 
only for the global energy budget, but to the energy budget of any climatic 
subdomains, with the crucial difference that the slowest internal time scale 
will be faster for, e.g., the atmosphere than for the ocean. 

Given the complexity of the processes inside the climate system and the 
fact that large cancelations are involved, the provision of a consistent or at 
least reasonable energy budget is in general far from being an obvious task, so 
that an extensive CMs intercomparison analysis is necessary to understand 
to what extent the actual CMs will conform to the case depicted in the first 
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panel of Fig. (1). 

In the case of renalyses, the bias in the diagnosed global energy budget 
is of the order of -10 Wm~ 2 , which is about one order of magnitude larger in 
absolute value (and wrong in the sign) than the reasonably expected imbal- 
ance due to the anthropogenic and natural forcings in the second half of the 
XX century. Surprisingly, the performance of these datasets for this specific 
benchmark does not improve even after 1979, which marks the beginning of 
the massive assimilation of the satellite data [69, 1]. 

Additionally, in spite of ever increasing and improving sources of obser- 
vative data from satellites on incoming and outgoing radiation at the top 
of the atmosphere (TOA), it is really hard to restrict the uncertainties on 
the observed yearly averaged energy balance of the planet to the order of 
1 Wm~ 2 [20, 63]. Very recently, it has emphasized how important is to 
correctly track, in models and observations, the Earth's energy for under- 
standing climate change and variability [61] 

Some studies have shown that numerical codes commonly adopted for 
the representation of the dynamics of geophysical fluids are not necessarily 
energy-conserving, as kinetic energy losses due to friction and diffusion are 
not exactly fed back in the system as thermal forcing, and that commonly 
adopted hyper diffusive schemes seem to be especially problematic in this 
regard [3, 31]. 

1.2 Transports 

When considering global energy balances the CS is reduced to a zero- dimensional 
system. We then make a step forward and tackle the problem of understand- 
ing how models represent the meridional enthalpy transport by the CS as a 
whole and by the atmosphere and the ocean separately. In order to achieve 
steady state, a net positive TOA energy budget in a given location must be 
compensated when long term averages are considered by a net divergence 
of the enthalpy flux. The TOA planetary budget has to a very good ap- 
proximation a zonal structure, with low latitude regions featuring a positive 
budget between the net input of solar radiation (controlled by incoming solar 
radiation and albedo) and the output of longwave radiation (controlled by 
surface temperature, water vapour concentration, and clouds). A negative 
budget is instead observed in high latitude regions [44]. Therefore, on the 
average, the fluid envelope of the CS transports large amounts of heat from 
low to high latitude regions in both hemispheres [53, 48]. See [11, 64] for 
recent theoretical treatments in a rather general setting. 

It is then found that in each hemisphere the maximum intensity of the 
meridional transport is in the long term equal to the absolute value of the 
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total energy imbalance in the low latitude region (or, equivalently, the energy 
deficit in the high latitude region), whereas its maximum is located at the 
transition latitude where the TOA energy balance vanishes [44, 48]. There- 
fore, the location and the intensity of the maximum of the transport provide 
fundamental information on the properties of the CS. 

Earlier reconstructions of the meridional transport using interpolations of 
radiosonde data - with unavoidable uncertainties due to the uneven spatial 
coverage and homogeneity issues of the datasets - suggested a comparable 
contribution to the meridional transport by the ocean and by the atmosphere 
[44]. Instead, more recent data suggest that, especially in the southern hemi- 
sphere, the atmospheric transport definitely plays a dominant role [58]. 

On the other hand, large scale organized atmospheric and oceanic motions 
are result of the mechanical work (then dissipated in a turbulent cascade) 
produced by the climate engine thorugh the Lorenz energy cycle, which is 
fuelled essentially by the presence of temperature differences across the CS 
[26]. Recently, a rather complete picture of the Lorenz energy cycle in a 
simple yet rigorous thermodynamic terms has been provided [19] and linked 
to estimates of the turbulent entropy production of the planet [29]. 

Many properties of the CS play an important role in determining the 
properties of the meridional atmospheric and oceanic transport, the main 
being the meridional gradient of the planetary albedo. Moreover, on Earth 
the mechanisms of transport are rather different in the northern and southern 
hemisphere. In the northern hemisphere, considerable contributions to the 
meridional transport are given, separately, by the transient and stationary 
waves of the atmosphere as well as by the ocean. Instead, in the southern 
hemisphere, most of the transport is performed by transient atmospheric 
waves only. Nonetheless, the maximum value and the location of the peak of 
the transport are very similar in the two hemispheres [58]. 

In fact, some basic features of the transport profile seem to be symmetric 
with respect to the equator as they are constrained by basic fluid dynamical 
laws and by the geometrical properties of the system. In a classical paper, 
Stone [54] showed that in a planet the climatological total meridional trans- 
port profile (and especially the position of its peak) is determined, to the first 
order of approximation, by spherical geometry and latitudinal dependence of 
TOA incoming radiation only. See also [9] for a recent re-examination of this 
brilliant idea. 

Since we still do not have a comprehensive rigorous theory able to express 
the peak location and strength of the meridional transport and to accomodate 
the variety of processes contributing to it, it seems crucial to test the degree of 
realism and consistency of CMs in representing these basic climatic features. 

The direct calculation of the transports from 3D climate fields requires 
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very serious diagnostic tools and high-resolution data (in time), since quadratic 
quantities are involved. Nevertheless, under the assumption of steady state 
conditions and of no biases in the energy balances of the whole CS, of the at- 
mosphere and of the ocean, it is possible to rigorously reconstruct the merid- 
ional enthalpy transport by atmosphere and ocean from the zonal averages 
of the heat fluxes at the upper and lower boundaries of the fluid envelope 
[44] . Even in the presence of biases in the energy budgets - if these are small 
enough with respect to the latitudinal variability - such a reconstruction tech- 
nique can be safely adopted with minor corrections, as discussed in [5, 58], 
where non steady state conditions like those of the present climate and rather 
heavily biased energy balances are analyzed. The meridional transport pro- 
files provide a very strong characterization of the CMs, and it seems crucial 
to intercompare their performances in order to understand how well they 
agree in the representation of such crucial climatic properties. Moreover, the 
analysis of the response to CO2 increase of the meridional transport probably 
provides the most basic way of investigating the adjustment of the dynamics 
of the fluid components of the CS. 

1.3 Goals and Structure of the Paper 

The overall goal of this paper is to present a complete review of the state- 
of-the-art of climate modelling in terms of their representation of the CS 
energetics and of the energetics of ocean and atmosphere seperately. The 
paper also wishes to stimulate lines of activities aimed at improving the 
reliability of CMs, which is a task of the uttermost urgency in the context 
of the preparation of the future IPCC reports. The paper is organized as 
follows. In Section 2, starting from the Eulerian energy balance equations, 
we briefly recapitulate the equations determining the energy budget of the 
CS and of its subdomains and show how the meridional transport can be 
inferred. In Section 3 we present the data considered in this work and discuss 
the methods used in their processing. In Section 4 we present our results, in 
Section 5 we discuss our findings, and in Section 6 we give our conclusions 
and provide perspectives for future work. 

2 Energy Budget in the Climate System 

In this section we review the derivation of the zonally and globally integrated 
energy balance of the CS and of its subdomains and show how, assuming 
steady state conditions, it is possible to infer the meridional enthalpy trans- 
port by the atmosphere, by the ocean, and by the system as a whole. We 
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also discuss how to derive such transport properties when energy budgets are 
not perfectly close, due to actual deviations from steady state or to spurious 
biases in the data. 

The evolution equation in Eulerian form for the total energy density e 
in the fluid component of the Earth system, given by the sum of the kinetic 
and moist static potential components, can be written as [44]: 

3 pc 

JL. = -V ■ (pec) - V ■ F R - V ■ F s - V ■ F L - V • (pc + r ■ c) (1) 

where p is the density, c is the velocity vector, Fr, Fs, and Fl are the 
vectors of the radiative, turbulent sensible and turbulent latent heat fluxes, 
respectively, p is the pressure and r is the stress tensor. It is useful to rewrite 
the previous equation as: 

^ = -V • J h - V • F R - V • F s - V • F L - V • (r • c) (2) 

where we have introduced the total enthalpy transport Jh = (pe+p)c = phc 
where the standard definition of enthalpy h = e + p/p is adopted [10, 23]. 
Note that the more commonly the expression energy transport is adopted 
[44]. As discussed in [2], Eq. (2) expresses a local balance corresponding to 
the first law of thermodynamics for an open system, and enthalpy (rather 
than energy) transports enter the picture because exchange of mass takes 
place. Note also that since all terms on the right hand side of Eqs. (1-2) are 
written as convergence of vector fields, by using Gauss' theorem it is easy 
to obtain an expression for the derivative of the total energy content of a 
subdomain of the CS in terms of the fluxes at the boundary. 

Expressing the divergence operator in spherical coordinates, and taking 
the usual thin shell approximation so that R + z f» R with R the radius of 
the planet, we obtain: 



dpe 1 dJh\ 1 d J hp cos <p dJj 



Itz 



+ 



dt Rcosp d\ R cos <p dip dz 

-§£ F R* + F sz + F Lz + (t-c) z ), (3) 

where in the last term of Eq. 3 we have retained only the vertical components 
of the radiative and turbulent fluxes, and of the contributions due to frictional 
stresses, since they are largely dominant by various orders of magnitude [44]. 

If we want to characterize the large-scale properties of the CS, we need 
to analyze the properties of the temporal and spatial averages of the energy 
budget. It is particularly relevant to exploit Eq. (3) to study the interaction 
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between the various climatic sub-systems. Let f2 be a subdomain of the CS. 
For our purposes, Q will refer to the atmosphere (Q=A), the ocean (Q=0), 
or the land subdomain (Q=L), as well as the complete system (Q=T). Given 
a generic field a=a(t, A, <p, z), we denote as mass weighted vertical integral of 
a over Q the integral 

Aa(t t X,<p)= f pa(t,X,<p,z') dz' (4) 



2(, 



where Zb and 2^ are respectively the lower and the upper boundaries that 
delimit Q. Note that where the thickness of the subdomain Q is vanishing - 
e.g., Q = O and the coordinates (A, tp) correspond to continental areas - Eq. 
(4) is still well-defined as A^it, X,ip) results automatically to be zero. We 
can then define the zonal integral: 

r-2-K 

[An](t,(p)= A n (t,X',ip)R cos ipdX' (5) 
Jo 

In order to characterize the 1-dimensional structure of the system, we inte- 
grate zonally and vertically Eq. 3. For the entire planet (Q = T) we have 

1 dT T 

[Et} = ~ 1i -^ + [F b ]toa (6) 
where [Fr\toa is the zonal mean of Fr evaluated at TOA and 

r2iT pTOA 

T T (t,<p)= / J hip {t,X , ,ip,z J )Rcosipdz , dX' (7) 

JO Jz b 

is the total transport of enthalpy realized through the latitudinal wall at 
latitude ip, being Zb the bottom of the ocean. Similarly, we can derive the 
equations for the atmosphere, ocean and land: 

[E A ] = + [Fr]toa + [F R + F s + F L ] surf>A ; (8) 

1 cfT 

[Eo] = -r-q£- ~ i F R + E s + F L ] surW - (9) 

[E L ] = -[F r + F s + F L ] surfyL] (10) 

where the subscripts in the last terms on the right hand side of Eqs. (8-10) 
refer to the boundary considered. Moreover: 

T A (ip) = / / J h(p (t,X',(p,z')R cos cpdz'dX' (11) 

JO Jz s 
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and 

To(<p)= [ f ' J hv (t,\',<p,z , )Rcoa<pdz'd\' (12) 

JO J z b 

are the atmospheric and oceanic total transport through the latitudinal wall 
at latitude tp, respectively. Note that, since land is not a fluid subdomain of 
the system, the large scale enthalpy transport through this medium is entirely 
negligible, so that no transport term is present in Eq. (10). Note also that 
we have neglected the transfer of mechanical energy from the atmosphere 
to the ocean since it is smaller by some orders of magnitude than the other 
energy exchange terms [44] . Such a term is instead of great importance when 
considering the driving mechanisms of the large scale ocean circulation [22]. 
Note that the fact that this term is so small compared to the thermal energy 
fluxes is related to the very low efficiency of the CS as a whole [29] . 

If we further integrate Eq. (6) and Eqs. (8-10) over latitude from — ir/2 
to 7r/2, we obtain the global energy budget and the energy budgets for the 
atmosphere, the ocean, and the land as follows: 

{E T } = {F R } TOA , (13) 

{Ea} = {Fr}tOA + {Fr\ SU t},A + {Fs}surf,A + {-£l}sw/,A> (14) 
{Eo} — —\Fn}surf,0 + {Fs}surf,0 + {Fl} surf,Oi 

(15) 

{El} = —{Fr] SU t},L + {Fs}surf,L + {FL}surf,L, (16) 

where we have defined: 

{X}(t)= /"^ \x\{t,<p)Rdip. (17) 

J-w/2 

for a generic variable [X] (£,</>). 
2.1 Time averages 

If the system is forced with the usual seasonal cycle and the internal and 
external parameters of the system are not time-dependent, we expect that 
after a suitable transient, steady state conditions are realized. As different 
domains of the climatic system have different internal time scales, the realiza- 
tion of steady state conditions occurs in a large variety of time horizons when 
the system is initialized from arbitrary initial conditions or when parameters 
are suddenly altered. In general, when time-dependent forcing in the form 
of parametric variations are present, a climatic subsystem is in quasi-steady 
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conditions when its internal time scales are much faster than the charac- 
teristic time scale of the parametric change. In this case, the subsystem is 
driven by the parametric modulation of the forcing. The atmosphere and 
land subdomains have a much shorter time scale than the ocean. The atmo- 
sphere and the land can be considered to be balanced over time scales of less 
than one year, while the ocean features time scales ranging from some years 
to several hundreds of years. The current climate radiative forcing (and, 
arguably, those projected under the SRESA1B scenario) are slow enough 
that the atmosphere can be expected to be in quasi-steady-state conditions, 
whereas the ocean lags behind [51]. Following [47], we can define the atmo- 
sphere and land subdomain as thin climatic subdomains, while the ocean is a 
thick subdomain. The presence of interacting subsystem with very different 
characteristic time scales is one of the key elements of complexity of the CS. 

Under steady state conditions, the long-term average of the derivative 
of any climatic fields needs to vanish, because otherwise trends would be 
present. When long term averages are applied to Eqs. (13-16), we obtain the 
general formula: 

{E n } = {Bn} = 0, (18) 

where the overbar indicates the time averaging operation, and Bq is the 
algebraic sum of the fluxes entering the domain thorough the upper and 
lower boundaries of the domain Q. Instead, the time-average of the zonally 
and vertically integrated energy equation for a generic domain Q can be 
written as: 

[EuM^-^ + [Bu]. (19) 

Under steady state conditions, the left hand side of Eq. (19) vanishes, so 
that the average transport can be computed as 

T a (<p) = - [B Q ](tp')Rd<p'. (20) 
J ip 

Since basically no substantial large scale enthalpy transport can take place 
through a solid medium, in stationary conditions the surface energy fluxes 
over land are expected to be locally vanishing when long term averages are 
considered. Instead, non-vanishing sustained local energy imbalances are 
possible at the atmospheric and oceanic interface and at TOA, since they 
are compensated by the divergence of the enthalpy transport realized by the 
motions of the system. 

Geometry imposes that we should have zero meridional transport at the 
poles. In agreement with that, Eq. 20 imposes Xh=0 at ip — tt/2 (North 
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Pole). Moreover, in steady state conditions the energy budget {Bq} is ex- 
pected to vanish, so that, consistently, Tq is zero also at tp = —ir/2 (South 
Pole). 

As suggested in [5, 58], the procedure of inferring meridional transports 
from zonal energy balances can be considered robust enough also in the 
case of non-steady state conditions and inhomogeneous information, as when 
gathering from various sources observational estimates on the energy balance 
in present climate conditions. The condition for the applicability of the 
inference is that latitudinal variability of [Bq] (</?) should be much larger than 
{Bq}/(2ttR), which is easily met in terrestrial conditions for the outputs of 
climate simulations (see below). Nevertheless, in order to have vanishing 
transport also at tp — —tt/2 we need to correct Eq. 20. The simplest way to 
accomplish that [5] is to redefine Tq as follows: 

5RvO = - f l2 \M orr ^')Rd^\ (21) 

where 

cow 

[B Q ] (ip) = [B Q }(p)-{B n }/27rR. (22) 

This amounts to assuming a constant flux correction throughout the surface. 
Different approaches, based, e.g., on weighting differently various latitudinal 
belts on the basis of the a priori knowledge of the data uncertainties, provide 
rather similar outputs when biases are overall small as in our case. Therefore, 
we will stick to the definition given in Eq. (21) throughout this paper. 

3 Data and Methods 

The PCMDI/CMIP3 (http://www-pcmdi.llnl.gov/) is a CMs intercom- 
parison project run by the Lawrence Livermore National Laboratories (USA) 
freely providing the access to a vast dataset of outputs of standardized sim- 
ulations performed by state-of-the-art CMs under various past, present and 
future climate scenarios. The PCMDI/CMIP3 includes data from over 20 
CMs and has been a key contribution to the IPCC4AR, since it has provided 
a large community of scientists the possibility of analyzing extensively and 
systematically the performances of many CMs on the same test-beds, and 
has greatly stimulated the development of metrics to be used for intercom- 
parison and validation studies. Presently, the new initiative PCMDI/CMIP5 
(http://cmip-pcmdi.llnl.gov/cmip5/) is being set up along similar lines 
to support the preparation of the fifth IPCC Assessment Report. Data refer- 
ring to this new generation of climate simulations should become available 
in the second half of 2011. 
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In this paper, we are interested in analyzing the behavior of the CMs in 
two standard IPCC4AR scenarios: 

• pre-industrial scenario, where the CS is integrated with fixed param- 
eters and in particular with fixed CO2 concentration set at 280 ppm 
until steady state conditions with no apparent trend in the climatic 
variables are observed; 

• SRESA1B scenario, characterized by a constant increase of CO2 from 
2000 to 2100, when 720 ppm concentration is obtained, and kept sub- 
sequently constant. 

It would have been interesting also to compare two different steady states 
(e.g., preindustrial vs. 2 x CO2 conditions) in order to characterize the 
generalized climate sensitivities of the system, but unfortunately the 2 x CO2 
control run data provided in the PCMDI/CMIP3 datasets refer to slab ocean 
integrations, and therefore are of limited use. 

In spite of the overall purpose of PCMDI/CMIP3 of providing complete 
and standardized datasets, we have found some deficiencies in data avail- 
ability, which lead to difficulties in performing a complete and systematic 
intercomparison analysis. This calls for further improvements for initiatives 
following the magnificent efforts by PCMDI/CMIP3. A typical problem has 
been the lack of one or more surface fields for a CM in a certain scenario, 
which has restricted the analysis to TOA properties only. In some cases, also 
TOA fields are lacking, so that the corresponding CM is simply excluded from 
the analysis for that scenario. Moreover, some CMs have corrupted data for 
certain fields or incorrect metadata accompanying the datasets. Therefore, 
some data have bene directly obtained from the modeling centers. 

Other inconsistencies are related to the temporal window of the data 
provided for different CMs for a given scenario. Regarding the PI scenario, 
the various CMs present time series of data of different lengths. In our case, 
we have chosen 100 years as the standard length of the time series to be 
analyzed. However, some CMs present less than 100 years of data. We 
have decided to perform the analysis also on these CMs, operating on the 
available data. Regarding the SRESA1B scenario, not all the CMs present 
data for the XXII and for the XXIII century: some CMs' data stop on 2200, 
others even on 2100 (before stabilization of CO2 has begun). Again, we 
have operated on all the available data after 2100. Finally, we have taken 
the radical choice of excluding the models with flux adjustment since their 
energetics is intrinsically biased - they are in contact with an infinite heat 
reservoir. An exception has been made for the INM-CM3.0 model since the 
flux adjustments are present only in a very limited oceanic region, so that no 
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substantial impacts are expected at global level and for zoanlly integrated 
quantities. An overview of the CMs considered in this study and of the data 
availability of the scenario runs is given in Table 1. Each CM is labelled by 
an identifying number: this convention is used throughout the text and in 
the figures. 

In the optimal scenario, we have considered for both climate scenarios 
the fields of short and longwave radiation fluxes at TOA and at the sur- 
face, plus the fields of latent and sensible heat fluxes at the surface. Ra- 
diative fluxes as well as latent and sensible heat fluxes are expressed in 
Wm~ 2 . Moreover, we have taken into account the Land-sea mask, which 
quantifies the portion of surface covered by sea and by land. Depend- 
ing on the considered CM, the values of the land-sea mask can be either 
boolean (0 = sea and 1 = land), or present continuous values between 
and 1, thus allowing for mixed coverage. Fields relative to the energy 
fluxes are 2+1 dimensional, being time series of bidimensional fields, while 
the land-sea mask is a fixed bidimensional field. Since we are interested 
in yearly averages, we have reduced the computational and storage bur- 
den of our facilities by downloading only pre-processed monthly mean fields, 
and then used them to compute the annual means. Data have been down- 
loaded in netcdf format, have then been imported in MATLAB 7.0® en- 
vironment with the freely available CSIRO netCDF/OPeNDAP software 
(http://www.marine.csiro.au/sw/matlab-netcdf.html), and have even- 
tually been processed with in-house scripts. 

Note that, in agreement with the actual estimates for the present climate 
[44], the surface average of the transfer of mechanical energy from the at- 
mosphere in the ocean is in all models and in both scenarios positive and 
of the order of 10~ 2 Wm~ 2 at most. This terms is several orders of mag- 
nitude smaller than the absolute value of the other terms contributing to 
the atmospheric and oceanic energy budgets. The tiny bias introduced by 
neglecting the transfer of mechanical energy does not, as we shall see below, 
add much to the analysis of the energy budgets, so that we can safely neglect 
this contribution. 

3.1 Statistical Estimators 

In the pre-industrial control runs, since the slowest external forcing to the 
system is the seasonal cycle (no variability of the incoming solar radiation 
is considered), steady-state conditions are expected to hold. Therefore, the 
time series of the yearly zonally integrated ([Bn])j and globally integrated 
energy balance {{B^})j (where j refers to the j th of the N years) are expected 
to have vanishing expectation values for each domain Q . Thus, for a given 
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Q, the yearly meridional transport one-dimensional fields (Tu(tp))j computed 
by applying Eq. (21) are expected to be statistically equivalent. 

For the XXII and XXIII century SRESA1B scenario outputs, the slowest 
forcing is also the seasonal cycle, so that generating the time series of yearly 
averaged values is a sensible operations. Nevertheless, as the system is read- 
justing to a new steady state, we may expect to see trends in the yearly data 
related to slow subsystems, i.e., the whole CS and the ocean subdomain. 
Actually, we find an approximate exponential behavior with decaying time 
which can be estimated as being of the order of 500-700 years, which nicely 
matches the characteristic time scale of the large scale ocean circulation [22]. 
For the time scales of our interest (order of 100 years), it can be treated 
safely clS db linear trend. 

The study of the global energy budgets relies on the analysis of the sta- 
tistical properties of the mean values and the interannual variabilities of the 
yearly time series ({i?^})^. The estimate of the yearly average of the energy 
budget is computed as the mean over the N years, and the estimate of the 
interannual variability is evaluated as the usual standard deviation of the 
time series. In the case the time series are not stationary, a linear trend has 
been removed before computing the standard deviation. 

Note that, since we are interested in the statistical properties of the cli- 
matology, we need to compute the range of uncertainty of the estimates. In 
each case, the statistical properties of the mean and of the interannual aver- 
age have been obtained by applying the standard block-bootstrap resampling 
technique [67, 34], where the resampling block has been chosen to be as long 
as the decorrelation time of the yearly time series. The 95% confidence inter- 
val has then been computed using a set of 1000 synthetic time series. Again, 
in the case of not stationary time series, a linear trend has been removed 
before applying the block bootstrap resampling. 

As expected, the lagged correlation of the yearly time series of the atmo- 
spheric and land energy budgets result to be statistically non significant as 
soon as the lag is larger than one year. Since the relaxation time of these 
subdomains is fast, no significant trend has been found also in the SRESA1B 
scenario data. Instead, the climatic and oceanic budgets time series feature 
a decorrelation time of the order of 5 years. Note that this is a simple way 
to see that the decadal climatic signals are mainly due to the ocean. 

In all cases, we find that the confidence intervals found with the block- 
bootstrap method agree very well with the asymptotic (for N ^> 1) results. 
The asymptotic theory says that in the case of decorrelated time series, the 
standard deviation of best unbiased estimate of the yearly average is approx- 
imately given by the interannual variability times the factor 2/y/N, while, in 
general, when a decorrelation time n is present, the factor is approximately 
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2y/n/y/N, which is still <C 1 if N ^> n (as in our cases). The standard devi- 
ation of the best unbiased estimate of the interannual variability is a factor 
of v2 smaller than that of the corresponding mean. 

The investigation of the meridional transport in the pre-industrial and 
SRESA1B scenarios is based on the analysis of the time series of the one- 
dimensional fields {Tn(ip))j. The best unbiased estimator of the yearly aver- 
age of the meridional transport is computed as the average of the N trans- 
port profiles obtained for each year, and the best unbiased estimate of the 
interannual variability has been evaluated as the usual standard deviation 
computed for each latitude. Again, trends have been removed when comput- 
ing the standard deviation. Since when analyzing the meridional transport 
we are not testing a zero hypothesis, as in the case of the global energy bud- 
gets, we are less interested in obtaining detailed confidence intervals of our 
best estimates. Moreover, since we are dealing with ID fields with obvious 
spatial correlation, the applicability of the above mentioned block bootstrap 
resampling technique is more problematic. Therefore, we will present conser- 
vative estimates of the uncertainties by using the standard deviation of the 
time series as proxy of the uncertainty of the best estimate of the mean. As 
we shall see below, the interannual variability results to be so small that no 
piece of information will be lost with this procedure. 

4 Diagnostics of the Climate Models 
4.1 Global Energy Budget 

Results are presented in terms of intensive quantities, so that the CS (Fig. 2) 
and atmosphere (Fig. 3) energy budgets are divided by the total surface of the 
planet, whereas the ocean (Fig. 4) and land (Fig. 5) budgets are divided by 
the total ocean and land surfaces, respectively. We wish to underline that the 
approach of studying surface energy fluxes is thermodynamically equivalent 
to but computationally much more efficient than analyzing the variations 
of the total energy of the geophysical fluids. The latter approach is quite 
unfeasible when an extensive CMs intercomparison analysis is performed, 
because it would require the analysis of three-dimensional time-dependent 
fields for each CM. 

We present for each domain the scatter plots of the 95% confidence inter- 
val of the mean value and of the standard deviation of the yearly averaged 
intensive energy balance in preindustrial conditions (first panel). We present 
also a scatter plot of the 95% confidence interval of the variations of the 
mean value and the standard deviation from preindustrial to climate change 
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conditions (second panel). The values includued in the first panel have been 
obtained with the block-bootstrap technique using 100-year long time se- 
ries, selecting as a general rule the last 100 years of the data available in 
the datasets. In the case of SRESA1B simulations, we have considered the 
2101-2200 data, and, when available, the data referring to the XXIII century. 
The 95% confidence intervals have been centered on the best unbiased sta- 
tistical estimators and their half-width computed as twice the corresponding 
standard deviation computed as described above. 

An immediate outlook on the problem of simulating physically sensible 
statistical properties of the energy budget is given in the second panel of Fig. 
(1), where we present the yearly time series (after application of a smoothing 
filter of 5 years for purely graphical reasons) of the TOA globally averaged 
energy budget for the GFDL2.1 model (CM 9). The preindustrial control 
run is depicted with the black line. The run performed under varying CO2 
concentration is depicted with the red line. Such run encompasses the XX 
century simulation (started in 1861 from year 41 of the pre-industrial run) 
and the XXI century portion of the SRESA1B simulations (started from 
the end of the XX century simulation), is depicted in red line. Finally, 
the blue line refers to the XXII and XXIII century SRESA1B run, where 
the CO2 concentration is held fixed at 720 ppm. The energy balance of 
the control unperturbed run is biased towards a positive value of about 0.5 
Wm~ 2 . After the energy budget increases as a result of the genuine heating 
due to the enhanced greenhouse effect due to increasing CO2 concentration, 
the energy budget decreases towards smaller values as the system relaxes 
towards the newly established steady state. The relaxation process occurs 
over a time scale of few centuries (not well resolved by the data) and can be 
safely attributed to the adjusting process of ocean state. Moreover, it is not 
clear whether the final energy budget bias will be analogous to the one the 
CM features in pre-industrial conditions. Therefore, as discussed above, it is 
in principle not trivial at all to estimate the actual net total heating of the 
planet due to the transition from 280 ppm to 720 ppm CO2 concentration. 

4.1.1 Pre-industrial scenario 

In preindustrial conditions, most CMs feature a positive TOA yearly averaged 
energy budget, and in not even one CM the confidence interval intersects zero 
(see Fig. 2, first panel). Therefore, we can say that all CMs have a spurious 
behavior, which is apparently not compatible with steady state conditions. 
This is not the result of an adjustment in process, since trends are wholly 
negligible, but is instead a real bias. Threfore none of the CMs seems to 
conform to the idealized behavior depicted in Fig. (1). The imbalance is 
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typically below 2 Wm~ 2 in absolute value, apart from 3 outliers (CMs 7, 
10, and 15) featuring positive imbalances ranging between 2.5 and 5 Wm~ 2 . 
Models 4, 12, 13, 14 and 16 give the best performances, as they feature 
biases smaller than about 0.2 Wm~ 2 in absolute value. In principle, an 
average positive inbalance of 1 Wm~ 2 at TOA, if no net flux across the land 
surface and the ocean bottom occurs, is consistent with a drift of the energy 
content of the fluid envelope of the planet of about 1.6 • 10 22 J /year. 

Very large discrepancies exist among CMs in the interannual variabil- 
ity of energy budget, since the computed standard deviations of the CMs' 
yearly energy budgets range between 0.1 and 0.6 Wm~ 2 , with typical values 
of about 0.3 Wm~ 2 . The presence of such a wide span of results requires 
some attention, as it suggests the existence of inter-model differences in the 
intensity and/or time scales of the negative feedbacks regulating the estab- 
lishment of an (approximate) energy balance over multiyear time scales. We 
find no obvious relation between the intensity of the interannual variability 
and the intensity of the bias of the energy budget. 

In order to disentangle the possible causes for global energy imbalances 
in each CM we need to analyze the energy budgets of the atmosphere, of 
the ocean and of land separately. Common sense would suggest under the 
hypothesis that the CM is energetically consistent, that TOA spurious energy 
imbalances are basically due to biases in the surface energy budget of an 
ocean still not close enough from steady state. Under this hypothesis, we 
should have virtually vanishing energy budgets of the atmosphere and of 
the land components, as they are thin components of the CS and they are 
quasi-steady state. 

What we instead find is that the biases in the energy budget of the at- 
mosphere - Fig. (3) - and land - 5 - are of the same order of magnitude of 
the bias of the energy budget of the ocean - Fig. (4). Therefore, it seems 
physically unreasonable to interpret the imbalances of the CS as a whole and 
of its subdomains as being only (or mainly) due to still persistent transient 
conditions. Note that an imbalance of 1 Wm~ 2 for the atmosphere corre- 
sponds to a staggering drift of about 3 K/year of its average temperature. 
Since this is obviously not observed in the data, it is clear that some funda- 
mental issues in the dynamics of the CMs need to be addressed. Therefore, 
it is clear that the observed biases in the global energy budget of the CMs 
cannot be explained as a byproduct of the fact that the ocean component 
has not reached a true steady state. 

An explanation for the energy imbalance of the atmosphere can be tracked 
to the lack of closure of the Lorenz energy cycle [25, 26]. In fact, CMs in ge- 
neeral do not have a detailed parameterization scheme for the reintroduction 
as thermal forcing of the energy dissipated by viscous dissipation, interaction 
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with the boundary layer, cloud processes, and diffusion. For a specific study 
in this direction please refer to [3], while an example of prototypical value 
can be found in [31]. If in the considered medium the kinetic energy is just 
lost by dissipatioon and not re-injected in the system as heat (as also in the 
simple case of Boussinesq approximation), we expect a compensating spuri- 
ous positive energy budget, which then does not translate into a temperature 
drift of the atmosphere, but actually compensates a ghost energy loss and 
keeps the system at steady state. 

Land modules typically feature positive energy budgets, whereas in steady 
state conditions the land component of CMs should see, when long term av- 
erages are considered, vanishing energy balance not only in the global mean, 
but in each grid point, since no compensating horizontal transport is possi- 
ble in a solid medium. Therefore, land modules must feature spurious water 
budget in the land climate modules [36], due, e.g., to relaxation schemes in 
the bottom boundary conditions. Otherwise, they must present biases in the 
treatment of phase transitions, like those related to the common procedure 
of calving ice and snow accumulating beyond given thresholds, as done, e.g., 
in [33]. In CM 22 The spurious energy sink associated to ice calving has been 
estimated to amount to about 0.3 Wm~ 2 when globally averaged [60]. 

We note that the energy balance of the atmosphere (and in principle also 
of the ocean where, actually, some contribution from a slow temperature drift 
cannot be ruled out) of the majority of CMs is positive, thus suggesting that 
typically there is an underestimation of the thermal energy produced at the 
end of the turbulent cascade dissipating the kinetic energy. Moreover, most 
land modules also feature a positive budget, suggesting that energy is lost 
because of inconsistencies in the treatment of phase transitions and water 
and heat fluxes. 

In all CMs the standard deviation of the annual energy budget time se- 
ries is much smaller for the land and atmosphere components than for the 
oceanic component, whose interannual variability is similar to that of the 
TO A budget. This confirms that the ocean acts as an integrator and plays 
a major role in controlling the variability of the global energy budget over 
multiannual time scales. 

The performance of the various CMs wildly differs from case to case, 
and, typically, for a given model, the quality of the representation of the 
energy budget of the various subdo mains is not consistent, thus pointing to 
the modular structure of the CMs. In any case, it should be noted that the 
biases we find for the atmosphere and ocean subdomains are at least one 
order of magnitude larger than the mechanical energy transfer term we have 
neglected, which implies that our approximation is consistent with the scope 
of our analysis. 
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We hereby highlight some specifically interesting features. We first note 
that CMs 4, and 16 feature a energy balance of the atmosphere and of the 
ocean accurately close to 0, with the land modules of introducing only a 
minor bias in the global energy budget. As opposed to that, CM 7 features 
a well-balanced land component, but both the atmosphere and the ocean 
feature strong positive biases, indicating that the energetics of the fluid media 
requires attention. On the other hand, the excellent performance on the 
global budget of CM 14 comes from non-negligible compensating errors in 
the sea, atmosphere and land modules. The poor performance of CM 15 in 
terms of global energy budget seems to derive from inconsistencies in the 
representation of the atmospheric and of the land energy balance. Energy 
inconsistencies in the atmosphere seem to be very relevant also for CM 1 and 
CM 10. 

4.1.2 SRESA1B scenario 

While for each CM it is crucial to assess the energy bias in the pre- industrial, 
steady state conditions, globally and in each subdomain, it is also important 
to understand how the energy budget chang result of a prolonged 

radiative forcing. The change in the energy budget between the pre-industrial 
and the SRESA1B scenario are presented in the second panels of Figs. 2-5. 
Values are shown for data referring to 2101-2200 (just numerical index) and, 
where available, also to 2201-2300 (same numerical index with lower index 
'a'). 

In the SRESA1B simulations most CMs feature an increase in the total 
energy budget, ranging from about 0.2 Wm~ 2 to 2.0 Wm~ 2 with a typical 
value of about 1.0 Wm~ 2 . This can be interpreted as the fact that the system 
is still gaining energy from outside as steady-state conditions have not yet 
been attained. The relaxation process in action is made more evident by 
the fact that the change in the energy bias for the XXIII century is 20-30% 
smaller than for the XXII century for all CMs where longer integrations are 
available (4, 8, 9, 12, 18, 20, and 22). Note also that the changes shown in 
the second panel of Fig. 2b, which should be taken as robust signature of 
climate change, are of the same order of magnitude or even larger than the 
biases found in steady state conditions (see first panel of Fig. 2a). 

Disentangling this dependence of the bias on the climate state among the 
various subdomains provides additional insights (second panels in Fig. 3-5). 
We find that the changes in the energy bias of the atmosphere are relatively 
small, ranging for all CMs between -0.1 Wm~ 2 and 0.3 Wm~ 2 . This is 
consistent with the fact that the atmosphere is in a quasi-steady state and 
supports the fact that the biases observed in pre-industrial conditions are 
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due to the internal dynamical processes and are relatively state- independent. 

Similarly, for almost all models the land energy balance changes by a little 
amount (less than 0.2 Wm~ 2 except CMs 15 and 16) between SRESA1B 
and pre-industrial conditions, again in agreement with the fact that we are 
dealing with a thin, fast system. Also in this case, the biases seem to depend 
weakly on the climate state. This provides further support to the fact that 
the discrepancies found in the pre-industrial CclSG £1X6, cLS for the atmosphere, 
related to structure of the land modules. 

Instead, the ocean acts as a heat reservoir for all CMs: the energy balance 
of the ocean is consistently increased for all the CMs between 0.5 and 2.5 
Wm~ 2 in the SRESA1B scenario with respect to the pre-industrial one, and 
it is the key player in determining the change in the global budget shown 
in Fig. 2. Note that values reported in Fig. 4b are larger than those given 
for the TOA in Fig. 2b because the ocean covers about 70% of the Earth 
surface: the extensive values of the change in the energy uptake are, instead, 
rather similar. The ocean is still relaxing to its steady state and is genuinely 
taking up heat. The relaxation to the steady state is confirmed by the fact 
that, as for the TOA budget, the values referring to the XXIII century are 
smaller than those referring to the earlier century by about 20-30%. 

The importance of discussing the energy budgets for each subdomain sep- 
arately can be emphasized by the following argument. By checking only the 
first panel of Fig. 2, we would conclude that there is a relatively large group 
of CMs (e.g., 4, 9, 11, 12, 13, 14, 16, and 17) featuring comparable perfor- 
mances in representing correctly the steady state properties of the CS. The 
analysis of the second panel of Fig. 2 does not allow us to discriminate among 
them, as we know that a non-zero response has to be expected. Instead, we 
observe that CMs 4 and 16 consistently provide among the smallest biases 
in energy budget of all subdomains in pre-industrial conditions and feature 
very small changes in the bias for the atmosphere and land components also 
in SRESA1B conditions, as theoretically expected since these are thin sub- 
systems. Therefore, we conclude that these CM are the best in terms of 
representing consistently the energy budget of the CS. Instaed, it is quite 
interesting to note that in the case of CM 9 a global energy balance close to 
results from an almost perfect cancellation o fthe atmospheric and oceanic 
biases. Unfortunately, we do not have enough data to assess the quality of 
CMs 11, 12 and 14 in all subdomains. 

4.2 Meridional Enthalpy Transport 

Composite plots of the profiles fi((T n ((p)) j) (Q = T, Q = A, Q = O) obtained 
for the CMs analysed in this study are shown in Fig. 6a for pre-industrial and 
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in Fig. 6b for SRESA1B conditions. Individual CMs results are not singled 
out (see below). All CMs feature a very low interannual variability on the 
meridional transport profiles, with absolute values of the local coefficients of 
variation (ratio of the standard deviation and of the mean) of the order of 
1% except near the region where meridional transport changes sign. There- 
fore, the variability is not reported in the figures. In each scenario, CMs 
agree on the most important features of the meridional transport, at least 
qualitatively, and no radical differences emerge between the two scenarios. 
This shows that the impact of the (robust) CO2 increase is to modulate the 
global circulation and not to disrupt it. We clearly see that the fluid envelope 
transport enthalpy from the low- to the high-latitude regions, with the effect 
of decreasing the meridional temperature gradient and of producing entropy 
[44, 43]. 

In both scenarios, for most CMs the atmospheric energy transport is, in 
agreement with today's climate [44, 58, 8], stronger in the southern hemi- 
sphere than in the northern hemisphere. In both hemispheres the peak of 
the atmospheric transport is located at the mid-latitudes, where we have the 
maximum of the baroclinic activity. Most CMs feature a slight southward 
atmospheric transport across the equator, probably due to the fact that the 
Hadley cell is not perfectly symmetric, since the Intertropical Convergence 
Zone is located north of the equator [45]. The atmospheric transport profile 
is very smooth, thus showing a strong dynamical link between the Hadley 
circulation and the mid-latitude atmospheric variability [62]. 

The ocean features a strongly asymmetric meridional transport, which 
mirrors the asymmetry in the oceanic circulation between the southern and 
the northern hemisphere [22]. In both hemispheres most CM feature two local 
maxima of transport, corresponding to the gyres. The larger one peaks in 
the low latitudes, while the secondary one peaks in mid-to-high latitudes. 
Overall, transport is stronger in northern hemisphere, with vanishing or 
weak northward cross-equatorial transport. The structure of the (weaker) 
meridional transport in the southern hemisphere is somewhat more com- 
plex, as some CMs feature latitudinal bands of oceanic equatorward counter- 
transport. As discussed in [60], the representation of ocean transport in the 
southern hemisphere is problematic aspect of most PCMDI/CMIP3 CMs. 

Finally, the total northward transport is almost hemispherically specular, 
featuring a very smooth shape, not exactly vanishing at the equator for 
most CMs. Since the peaks of the atmospheric and of the oceanic transport 
are shifted by more than 20°, and the atmospheric transport is typically 
stronger by a factor of ~ 3 — 5, the peak of the total transport is slightly 
shifted equatorward with respect to the atmospheric one. The most relevant 
qualitative feature of the meridional transport profiles presented in Figs. 6a- 
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b match remarkably well with those inferred from observations of the present 
energy budget at TOA and at surface [58], thus implying that they constitute 
robust structural properties of the CS which are relatively well captured by 
the CMs. In particular, it is quite confuting that, in agreement with Stone's 
results [54], for all CMs the total transport profile is remarkably similar in 
the two hemispheres in spite of wide difference in the contributions provided 
by trasnient and stationary waves in the atmosphere and by the oceanic 
transport. 

From a quantitative point of view, the disagreements among CMs are 
large. In the pre-industrial simulations, the northern hemisphere total and 
atmospheric transport features discrepancies among CMs up to 10%, and 
up to 20% in the southern hemisphere. Agreement is worse for the oceanic 
transport, discrepancies being up to 40-50% in the northern hemisphere. 
In the southern hemisphere, where the intensity of the oceanic transport is 
lower, we have discrepancies up to 100%, and also the qualitative behaviour is 
less consistently represented [60]. In the SRESA1B simulations the situation 
is rather similar, with the notable difference that, even if the number of 
models included in the analysis is smaller, the models discrepancies on the 
ocean meridional transport are larger, probably as a result of the fact that 
transient time response of the ocean to CO2 is rather uncertain. 

Efficient metrics for the evaluation of CMs require the definition of simple 
but comprehensive indicators able to characterize the main features of the 
process under investigation. As an example, in a previous work we scaled 
down the Hayashi spectral density of mid-latitude waves to some integrated 
variables describing the overall spectral weight of qualitatively different waves 
[35]. In this case, we define a procedure to extract from the ID fields of 
meridional transports the essential information. Therefore, from each yearly 
averaged transport profile (Th(y?))j we extract four variables, namely the 
value of the peak of the northward transport in the northern hemisphere 
(T^ H )j, its latitude (<fn H )j, the peak of the southward transport in the 
southern hemisphere (Tn H )j, and its latitude {^ H )j- The latitudinal band 
between (fQ H )j and (<fQ H )j defines the region within Q from which a net ex- 
port of enthalpy is realized (except the case of southern hemisphere oceanic 
equatorward counter-transport), whereas (T^ H )j {(T^ H )j) describes the to- 
tal imbalance in the domain Q north of (south of) ( l fn H )j {(<Pn H )j)- 

The various CMs feature rather different horizontal grids, in terms of 
resolution and structure (e.g., regular lat-lon vs. Gaussian grids), so that 
the comparison of local indicators like those introduced above requires the 
remapping of all profiles (Tq (if)) j into the same latitudinal grid. Therefore, 
we have selected the regular 2.5° spacing from —90°N to 90° N as a common 
framework , since it approximately corresponds to the average resolution of 
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the CMs, and performed the remapping for all CMs using a linear interpola- 
tor. Results are robust with respect to the interpolating algorithm. 

4.2.1 Pre-industrial scenario 

We begin our analysis by focusing on the peaks of the total meridional trans- 
port in pre-industrial conditions for the northern and southern hemisphere, 
depicted in Figs. 7a-b, respectively. We first observe that in both hemi- 
spheres most CMs feature transports peaking at 5 to 6 Petawatts (PW, 1 
PW = 10 15 W), so with an effective range of about 20%, and with max- 
ima situated between 34° and 38°. On the average, the total transport in 
the sourthern hemisphere is typically slightly smaller by about 0.3 PW, in 
agreement with what found in [60] for present climate conditions on a sim- 
ilar set of CMs. The model-wise disagreements on the localization of the 
transport maxima are more serious in the northern hemisphere, while dis- 
agreements on the peak value of the transport are somewhat wider in the 
southern hemisphere. 

For each CM, the interannual variability of the position and the strength 
of the transport maximum is very small, so that typically each year the merid- 
ional transport peaks at the same grid point, and the transport maximum, as 
also discussed in the context of Fig. 6a, has weak year-to-year variations. This 
is hardly unexpected, since we are looking into one of the most robust fea- 
tures of the CS, which should be affected only by large climate fluctuations. 
We may guess that when the CS is close to tipping points, the enhanced, 
large scale, organized fluctuations might heavily impact also these indica- 
tors. An interesting result is that CMs tend to behave consistently in the 
two hemispheres. Typically, relatively strong (weak) transports are observed 
in both hemispheres, as in the case of CMs 10, 15, 16, 20, 22 (CMs 4, 5, 11, 
12, 17). Similarly, relative poleward (equatorward) anomalies in the position 
of the transport maxima are often consistent in the two hemispheres, as in 
the case of CM 17 (CMs 5 and 6). 

The agreement among CMs in the representation of the atmospheric en- 
thalpy transport in the northern (see Fig. 8a) and southern hemisphere (see 
Fig. 8b is also problematic. In both hemispheres the agreement among CMs 
is rather good on the position of the maxima (within one grid point), with 
a slight more poleward position of the maxima in the northern than in the 
southern hemisphere. Since we are dealing with latitudes of about 40°, we 
deduce that a great contribution to the transport is performed by baroclinic 
activity, and that the position of the peak coincides with that of the storm 
track. As in the case of the total transport, the interannual variability of the 
position and the intensity of the peak of the atmospheric transport is very 
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small for all CMs. 

Instead, quantitative differences in the intensity of the transport are quite 
pronounced in both hemispheres. In the northern hemisphere, most values 
are clustered within 1 PW around 4.7 PW, which means a spread of about 
20%. Such inconsistencies in the representation of the large scale properties 
of the atmosphere in the northern hemisphere agree with what discussed in 
[35]. In this paper, albeit referred to the present climate, it is shown that 
large discrepancies among PCMDI/CMIP3 CMs exist in the representation of 
the northern hemisphere winter atmospheric variability in the mid-latitudes, 
which is when and where the meridional energy transport is particularly 
effective. 

In the southern hemisphere, even if the total meridional enthalpy trans- 
port is typically lower, for most CMs the values of the maximum of the 
atmospheric transport tend to be somewhat larger than in the northern hemi- 
sphere, by up to 0.5 PW. This is in agreement with what found in [7, 8] for 
present climate data, and it is basically due by the fact that the baroclinic ac- 
tivity is stronger in the southern hemisphere [57]. The intensity span among 
CMs is more pronounced than in the northern hemisphere, reaching a value 
of 1.5 PW (about 30%), thus suggesting more uncertainty among CMs in 
the representation of the the meridional temperature gradient in the southern 
hemisphere. We also observe that some CMs feature extremely weak (CMs 
1, and 4) and extremely strong (CMs 15, 16, and 20) atmospheric transports 
in both hemispheres, which hints at the presence of some peculiar structural 
properties of the atmospheric circulation. 

When looking at the CMs performance in terms of representation of the 
ocean meridional enthalpy transport, the picture is rather different than in 
the two previous cases. We find results supporting the idea the the dis- 
crepancies among CMs in the representation of the ocean circulation - and 
especially that of the southern ocean - are much larger than in the case of 
atmospheric circulation. 

In preindustrial conditions, in the northern hemisphere the intensity of 
the maxima of ocean transport are between 1.3 PW and 2 PW, so the span is 
about 40% of the typical value. The positions of the peaks are not correlated 
to their intensity and are distributed between 15°iV and 25°iV, so with a span 
of 4 grid points. The maximum for CM 10 is relevantly shifted northward, 
up to around 30°iV. For each CM, interannual variability is small for both 
intensity and position of the maximum, thus supporting the idea that the 
peak of the oceanic transport in the northern hemisphere provides a robust 
characterization of the state of the oceanic system. 

As already seen in Fig. 6a, in the southern hemisphere the definition of 
the peak of the oceanic transport is more problematic due to the presence of 
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two well defined local maxima. We can see that for most CMs the peak of 
the oceanic transport is located in a relatively narrow band between 12.5°S 
and lO ^ (1 grid point span), while for CM 6 the peak is shifted poleward 
between 15° S and 20°5. This CM typically features a very weak transport, 
peaking typically around 10° S, but in few specific years the maximum of the 
transport is realized around 55° S. Therefore, the mean peak position has an 
intermediate value. Finally, the bimodality in the location of the maximum is 
very pronounced in CMs 1, 4, and 10, for which the (rather weak) meridional 
transport peaks in high latitudes quite often. A more detailed analysis of the 
distribution of the maxima accounting for the bimodality of the pdf of the 
locations would surely be of interest. Restricting ourselves to the CMs with 
a clearly located maximum, we can see that the span of the intensities is very 
large. We observe values ranging from 0.4 PW to 1.6 PW, so that the span 
is about 100% of the typical value. Instead, for each CM, the interannual 
variability on the value and position of the peak is rather small. The existence 
of such discrepancies among the observed climatologies is in agreement with 
findings of previous studies, which showed that the PCMDI/CMIP3 models 
feature marked difficulties in the representation of large scale energy processes 
in the southern ocean. 

4.2.2 SRESA1B scenario 

Rather interesting differences emerge when comparing the simulations per- 
formed with pre-industrial scenarios and those performed with increased CO2 
concentration. The first non-trivial information is that for each CM the po- 
sition of the peak of the total meridional transport does not change signifi- 
cantly, in both the southern and the northern hemispheres. The shift is at 
most of one grid point and, model-wise, no systematic (poleward or equa- 
torward) variation is found (not shown). The behavior of each CM is in 
agreement with Stone's results mentioned above [54] on the stability of the 
position of the peak of the transport independently, e.g., of the specific at- 
mospheric composition. Therefore, the latitudinal boundary separating the 
areas of the planet differing on the sign of the difference between the incom- 
ing and outgoing radiation is rather stable with respect to changes in the 
CO2 concentration. 

Changes in the intensity of the meridional transport are more interesting. 
In Fig. 10a we plot the change in the value of the peak of the meridional 
transport in the southern hemisphere vs. the change observed in the north- 
ern hemisphere. Note that in all cases the 95% confidence interval is very 
narrow. We have that most CMs feature positive variations in both hemi- 
spheres up to about 10% of the pre-industrial value (typical changes range 
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around 5%), which implies that in a warmer climate energy is redistributed 
more effectively. This contributes to the mechanism of polar amplification 
[51], as its additional heat is fed into the high-latitude regions of the globe. 
Given the smaller thermal inertia of the system, the effect in terms of sur- 
face temperature is expected to be stronger in the northern hemisphere. In 
CM 22 significant deviations from this behavior are observed because, while 
the change in the meridional transport in the northern hemisphere is almost 
negligible. This CM foresees a relevant reduction (of the order of 6%) in the 
poleward meridional transport in the southern hemisphere. Therefore, in this 
CM the climate as a whole becomes less efficient in warming up the Antarc- 
tic region, and an increase in the relevant meridional temperature gradient is 
expected unless compensating changes in the meridional albedo gradient oc- 
cur. Instead, CMs 7, 14 and 15 foresee a negligible change in the meridional 
transport in both hemisphere, thus featuring a very stable climate. 

When looking at changes in the meridional atmospheric transport, we 
again find that in both hemispheres (see Fig. lla-b) changes in the position 
of the peak are small for all CMs but with few exception systematically pole- 
ward. As the atmospheric transport peaks are related to baroclinic activity 
and storm track, this results seems to be (at least partially) in agreement 
with what found in [70], who observed in XXI century climate simulations a 
poleward shift of the core of the mid-latitude atmospheric variability, with 
a stronger signal in the southern hemisphere. Nonetheless, recent analyses 
have attributed the shift of the peak in the southern hemisphere due to spu- 
rious trends in the clouds cover [60]. The impact of CO2 increase is instead 
relevant when the values of the peak of the atmospheric transport is con- 
sidered. In Fig. 10b we show that all CMs feature consistent increases (up 
to about 12%) of the atmospheric transport in both the northern and the 
southern (where changes are typically larger) hemisphere, with the exception 
of CMs 14 and 15, which foresee only a negligible increase of the transport for 
the northern hemisphere. Our results agree with what shown in [17], where 
aggregated results for the CMs' ensemble mean are presented. 

The response of the ocean to CO2 increase is also of great interest. We 
observe that typical changes of the order of 0.1-0.2 PW (Fig. 10c). These 
constitute large relative changes in the ocean transport, ranging from 10% to 
over 50%. Half of the CMs foresee a decrease in the transport in the south- 
ern hemisphere, while a decrease in the transport is foreseen for all CMs in 
the northern hemisphere. The observed (partial) compensation between the 
changes in the oceanic and atmospheric transports is in overall agreement 
with the aggregated results presented in [17]. There is no obvious link be- 
tween the observed changes and the unperturbed values. Instead, in both 
hemispheres most CMs feature (typically small) equatorward shifts in the po- 
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sition of the maxima, typically of the order of half of a grid point. The largest 
equatorward shifts are realized for CM 18 in the northern hemisphere and 
CM 4 in the southern hemisphere, which featured among the most poleward 
position of the maxima in the pre-industrial conditions. 

4.3 Atmosphere-Ocean Transport compensation 

The Bjerknes compensation mechanism suggests a compensating role of the 
atmosphere and of the ocean in transporting heat poleward [4]. The mecha- 
nism has been the subject of detailed studies performed also on PCMDI/CMIP3 
CMs [49, 55]. In the previous subsection we have seen that for higher CO2 
concentration for basically all CMs the increase in the atmospheric transport 
is moderated by the decrease in the oceanic transport in both hemispheres. 
We test here whether a hemispheric compensating mechanism is in place 
also on smaller time scales by computing for both hemispheres the linear 
correlation between the yearly time series of the maxima of the meridional 
transport for the atmosphere and for the ocean. Results are shown in Tab. 
2. Note that the 95% confidence interval of the null hypothesis corresponds 
to that of an uncorrelated time series since the atmospheric transport has no 
memory. We find that no CM, in no scenario, in neither hemisphere features 
a statistically significant positive correlation between the two time series. In- 
stead, we find that the majority of the CMs feature a statistically significant 
negative correlation between the time series in both hemispheres and in both 
scenarios. CMs 4, 6, 7, 10, 16, and 20 seem to have an especially strong and 
consistent compensation mechanism. We see that in both hemispheres the 
negative correlation is often reinforced in the SRESA1B scenario, suggesting 
that the mechanism may become stronger in a warmer climate. This may be 
motivated by the stronger coupling due to intensified (thanks to Clausius- 
Clapeyron relation) latent heat fluxes under global warming conditions. The 
fact that the Bjerknes compensation mechanism is active in both hemispheres 
and is evident also when yearly time series are considered further suggests 
that it acts a local or regional scale. In fact, if the compensation were re- 
lated to large scale ocean-atmosphere coupled processes, the behavior would 
be totally different in the two hemispheres, given the strong asymmetry of 
the oceanic circulation. 

5 Discussion 

In this work we have tackled the problem of assessing the consistency of a 
large set of state-of-the-art CMs in their representation of the energy balance 
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of the whole CS as well as of its main subdomains (atmosphere, ocean, land). 
Moreover, we have discussed the statistical properties of their meridional 
enthalpy transports due to the fluid components. 

Our analysis has consider the preindustrial and SRESA1B integrations 
performed following the standards provided by the PCMDI/CMIP3 project, 
which has greatly contributed to the redaction of the IPCC4AR [51]. The 
preindustrial scenario has played the role of a control run with fixed CO2 
concentration and has been instrumental for defining a baseline for dealing 
with transient-free climate behavior and for setting the objects of our investi- 
gation in the well-defined context of the non equilibrium steady state systems 
[12]. The other scenario has been selected in order to understand how the 
energetics of the CS chang result of changes in the CO2 concentration 

We find that in the pre-industrial simulations no CM features a global 
energy balance rigorously consistent with the condition of stationary state, 
so that a driftless state is realized even if no global balance is realized. The 
performances of the various state-of-the art CMs are wildly different, both 
in terms of global budgets, and when budgets of the main subdomains are 
considered. When global balances are considered, CMs biases are positive 
in all cases, with values spanning between 0.2Wm~ 2 and 2Wm~ 2 , with a 
few CMs featuring imbalances larger than 3Wm~ 2 . Typically, this behavior 
is interpreted as resulting almost entirely from biases in the surface energy 
budget of an energy- absorbing ocean, which has not attained yet its steady 
state, with an ensuing almost unnoticeable spurious but harmless drift in the 
climatogical temperature thanks to the ocean's large heat capacity. Instead, 
the energy imbalances in the land and atmospheric subdomain are of the same 
order of magnitude (about 1 Wm~ 2 ) of those of the ocean, which rules out 
the usual interpretation. Instead, inconsistencies are mostly related to basic 
issues in the land and fluid modules of the CMs. Note that for some of the 
CMs good global balances result from compensating errors in the subdomains. 
Instead, CMs 4 and 16 perform consistently well in all subdomains, which 
suggest that their representation of global energy processes is more physically 
consistent. 

Most CMs feature spurious positive energy budgets for the land subdo- 
main, as a result of unphysical biases in water fluxes [36] and in the tretament 
of phase transitions. As an example, long-term accumulation or loss of 15 
mm/year of liquid water (which amounts to about 1.5% of the intensity hy- 
drological cycle) or 12 mm/year of ice, e.g., due to spurious water infiltration 
or ice calving, correspond to an energy bias of about 1 Wm~ 2 in absolute 
value. 

CMs feature biases in the atmospheric energy budget, which are typically 
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positive and with of the order of 1 Wm~ 2 . It this were related to a genuine 
energy input to the atmospheric subsystem, it would correspond to a stag- 
gering average temperature trend of 3 K/year. As such a trend is definitely 
not observed, the spurious energy budget must compensate an unphysical en- 
ergy sink (if the budget if positive) or source (if the budget is negative). The 
inconsistency in the energy budget of the atmosphere is mainly related to 
the fact that kinetic energy dissipated by various processes, including viscos- 
ity and diffusion, cloud parameterisation, and interaction with the boundary 
layer, is not exactly re-injected in the system as thermal energy. This has 
been analyzed in detail in [3], whereas an example of prototypical value has 
been provided in [31]. As an example, commonly used hyperdiffusion schemes 
feature mathematical properties that make them inadequate in terms of rep- 
resenting correctly the energy budget and introduce an artifical energy sink. 

Note that these inconsistencies may have relevant impacts in the global 
properties of the atmospheric and and of the oceanic circulation, as they 
alter the input-output terms of the Lorenz energy cycle. Recently, in an 
intermediate complexity CM it has been shown that, when the energetic in- 
consistencies of the fluid and of the land modules are fixed, the representation 
of the global energy balance in greatly improved, with an obtained reduction 
of the the spurious imbalance by about one order of magnitude [37, 38]. 

We can link the presence of a bias in the global energy balance of the 
system AE to a bias ATe in the radiative temperature Te of the planet. 
Using a linear expansion to the black body relation, we find that ATe / AE = 
— 1/(4<tT|) (« — 0.27K/(Wm~ 2 ) with usual Earth values). Moreover, if we 
consider the simple but effective parameterisation for the outgoing longwave 
radiation expressed as the linear relationship A + BT$ with respect to the 
surface temperature T$, we find that an energy bias AE is related to a 
surface temperature bias ATs/AE = —1/B, which, using the - typical - 
Nakamura et al. values [41], results to be « — O.QK/(Wm~ 2 ). Therefore, we 
have that a positive (negative) imbalance - which actually cures an energy 
leakage (source) in steady state conditions - is related to a spurious cooling 
(warming) of the system, both in terms of its thermodynamic temperature 
and of its surface temperature. As basically models tend to have positive 
energy budgets in pre-industrial as well as in increased CO2 conditions, we 
expect them to feature a cold bias. We believe this is an alternative way of 
rationalizing the long-standing problem of the cold bias of CMs attributed, 
instead, by Johnson [18, 19] to spurious entropy generation due to numerical 
noise under the condition of vanishing energy balances. 

Another surprising result worth mentioning is that the interannual vari- 
ability of the energy budgets is model- wise very different, with values span- 
ning almost and order of magnitude, both in the global case and when sub- 
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domains are considered. This calls for a detailed re-examination of the differ- 
ences between the strength of the feedbacks controlling the energy budgets 
in the various CMs. 

When looking at energy biases obtained with the SRESA1B scenario runs, 
we find that the atmosphere and the land biases change by a rather small 
amount (typically up to 0.2 Wm~ 2 ) with respect to the pre-industrial runs. 
Since we are dealing with thin subsystems with small heat capacity, this 
suggests that the energy bias in not heavily climate-dependent but actually 
related to deficiencies in structure of the module. Instead, the ocean energy 
budget increases for all CMs, in agreement with the fact that in transient 
global warming conditions the water masses genuinely take up heat since the 
system is not yet at steady state. It must be observed that the spread of 
the heat update is very large, going from about 0.2 to over 2 Wm~ 2 . Such 
large discrepancies contribute to explaining the relatively large disagreement 
in the simulated transient climate change at surface by the various CMs. 

The analysis of the inferred meridional enthalpy transports has instead 
provided us with a picture of how CMs represent the process of re-equilibration 
by which the fluid component of the CS decrease the temperature difference 
between the low and high latitude regions in both hemispheres. The quali- 
tative agreement amongst CMs is rather good, and the obtained meridional 
transport profiles resemble recent estimates obtained with observative data 
[58] and CMs diagnostics studies [60], thus confirming a fair representation 
of robust structural properties of the CS. 

It is important to note that the indirect calculation of the transport is 
a linear operation, since it involves the integration of the vertical fluxes of 
energy at the boundary of the fluid envelopes and at the TOA. Therefore, 
we need to consider only 2D diagnostic fields with the temporal resolution of 
the climatology we wish to compute (yearly data, in this case). As opposed 
to that, the direct calculation of the transports requires 3D diagnostics fields 
with high temporal resolution, since quadratic quantities are involved. There- 
fore, our approach requires storage and processing costs smaller by orders of 
magnitude than what needed with a direct calculation. While the indirect 
calculation is rigorously valid only if steady state conditions are precisely 
obeyed and the integrated energy balances are vanishing, the bias found in 
the CMs are small enough that only negligible errors are introduced. 

Since we are dealing with ID fields, the necessity to deal with simpler 
metrics has driven us to consider more manageable indicators. We have then 
analyzed the statistical properties of the intensity and position of the hemi- 
spheric peaks of the atmospheric, of the oceanic and of the total transport. 

In pre-industrial conditions, we have found that the CMs do not agree 
precisely on the estimate of the total transport driven by the CS from the low 
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to the high latitudes regions, as discrepancies of the order of 15-20% around a 
typical value of about 5.5 PW (5.2 PW) are found in the northern (southern) 
hemispheres. Recent studies have suggested that the PCMDI/CMIP3 CMs 
might slightly underestimate the total transport in the southern hemisphere 
because too little radiation is absorbed in the tropical regions, due to to a 
cloud bias [60]. When considering the atmospheric transport, we find that in 
the northern hemisphere the CMs feature discrepancies in the peak intensity 
of 15-20% around a typical value of about 4.7 PW, whereas in the southern 
hemisphere the peaks of the transport are typically larger by about 0.5 PW 
and the discrepancy among CMs is more evident. Model 4 features an atmo- 
spheric transport much weaker than all other CMs. The agreement among 
CMs is much better for the atmospheric transport when looking at the posi- 
tion of the peak, with disagreements of only one grid point. The atmospheric 
transport peak is situated between 40° N and 43° N in the northern hemi- 
sphere and 37° S and 40° S in the southern hemisphere. The peak of the total 
transport is shifted equatorward in the northern (southern) hemisphere by 
about 5° (3°) and the CMs spread is somewhat larger. This confirms that, in 
spite of the serious discrepancy in the mechanisms contributing to the merid- 
ional transport in the two hemispheres, the overall properties are very closely 
symmetric with respect to the equator. This is in excellent agreement with 
Stone's theoretical results [54]. It is important to note that for each CM the 
interannual variability of the position and intensity of the transport maxima 
is very small, so that we are dealing with a very robust climatic feature. 

The inter-model uncertainties on the position of the peak of the total 
transport are mainly due to the large disagreement among CMs on the posi- 
tion of the peak of the ocean transport and on its intensity. Previous analyses 
have shown that the ability of PCMDI/CMIP3 CMs in describing the ocean 
transport in the southern hemisphere is especially problematic [60]. All CMs 
agree on the fact that the oceanic transport is smaller than the atmospheric 
transport in both hemispheres, and that the oceanic transport in the southern 
hemisphere is much weaker that that in the northern hemisphere. Neverthe- 
less, large quantitative inconsistencies emerge: disagreement on its intensity 
is of the order of 50% in the northern hemisphere (with values ranging from 
about 1.3 to about 2.1 PW) and of the order of 100% in the southern hemi- 
sphere (with values ranging from about 0.4 to about 1.6 PW). All CMs also 
agree on that the peak of the oceanic transport is shifted equatorward with 
respect to the peak of the atmospheric transport in both hemispheres. In- 
stead, large quantitative disagreements exist on the actual positions of the 
peaks. In the northern hemisphere, for all CMs except CM 10, whose peak is 
shifted poleward, the ocean transport peaks between 15° N and 25° N. In the 
southern hemisphere, for most CMs the peak is typically situated between 
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10° S and 20° S, whereas CMs 1, 4, 6, and 10 feature a much more poleward 
shifted average position. As opposed to the case of CM 10 in the northern 
hemisphere, this is not related to a permanent feature of the oceanic circu- 
lation. In the southern hemisphere for all CMs the ocean transport has two 
distinct peaks, one located in subtropical and the other one in high latitude 
regions. In the case of CMs 1, 4, 6, and 10, the high-latitude peak is typically 
about as strong as the subtropical one, so that the distribution of location 
of the yearly maxima of transport is basically bimodal. 

When looking at the SRESA1B scenario, we find that all models foresee 
an increase - ranging from 2% to 12% - in the atmospheric heat transport 
in both hemispheres. It is apparent that the enhancement of the latent heat 
fluxes due to increased mean temperatures dominates the foreseen reduction 
of the baroclinicity. A slight poleward shift is found for the peak of the 
atmospheric transport in both hemispheres, in broad agreement with what 
found in [70]. We have to remark that recent analyses attribute the poleward 
shift in the southern hemisphere to spurious representation of the changes in 
the cloud cover in altered climate conditions [60]. 

Following [17], the increase in the intensity of the atmospheric transport 
may be interpreted as the result of a strong cancellation between the large 
increase in the poleward latent heat transport, due to the strong enhancement 
of the hydrological cycle, and the reduction of the sensible heat poleward 
transport, due to the decrease of the baroclinicity of the system. A recent 
publication [42] has shown rather convincingly using a simplified yet Earth- 
like CM that this compensation mechanism is in place on rather wide range 
of climate scenarios. 

When considering the ocean component's response to CO2 increase, we 
find that all CMs foresee a substantial decrease in the intensity of the ocean 
transport in the northern hemisphere (disagreement exists in the southern 
hemisphere), and, typically, a relatively small equatorward shift in the po- 
sition of the peak in both hemisphere. These changes partially moderate 
the observed increase in the atmospheric transport, so that the total climate 
transport follow closely but not exactly those of the atmospheric transport: 
the peak is not shifted systematically in either hemisphere (again, in agree- 
ment with Stone's argument) and the transport is typically increased for 
most CMs up to 10% in both hemispheres. 

The contrasting response to CO2 concentration increase in the atmo- 
spheric and oceanic transports suggests that the Bjerknes compensating 
mechanism [4] is in place on climatological time scales. The yearly time 
series of the maxima of the atmospheric and oceanic transports have been 
used also to test - in a very unsophisticated way - whether the Bjerknes 
mechanism is active also on shorter time scale. The answer is partially pos- 
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itive, as both in the northern and in the southern hemisphere several CMs 
feature a statistically significant negative correlation between the time series, 
the rest features statistically non-significant correlations, while no cases of 
significant positive correlations are detected. The mechanism seems to be- 
come stronger in warmer climate conditions, possibly as a result of stronger 
atmosphere-ocean coupling due to more intense latent heat fluxes. 

6 Conclusions 

These findings suggest that, in spite of the great successes of climate model- 
ing, serious efforts seem to be needed in terms of basic science for improving 
substantially CMs. This is of great relevance when devising strategies for 
the future direction to be taken in the development of numerical simulations, 
and is presently of great interest in the context of the preparation of the 5 th 
Assessment Report of the IPCC. While earth system science requires models 
of increasing degree of complexity able to simulate correctly the chemical and 
biological cycles of our planet and their changes in a changing climate, we 
cannot expect to improve our understanding and ability to describe the CS 
just by adding and coupling indefinitely new modules. Nor the advocated 
quantum leap in climate simulations can be obtained just by improving the 
computational capabilities [50]. In particular, the definition of strategies for 
improving the closure of the energy cycle in the fluid components of the CS [3] 
seems to be a crucial step, as important as taking care of representing consis- 
tently phase transitions and heat fluxes over land. Some CMs already adopt 
the strategy of re-injecting in the system the kinetic energy lost to dissipation 
as a uniform heating term. While this helps in reducing the energy imbal- 
ance, it creates regional temperature biases, since climate regions featuring 
almost inviscid, horizontal dynamics are overheated (e.g., the stratosphere), 
and defines a spurious source of entropy production, as heat is moved down 
the gradient of the temperature field [19, 29]. 

It is apparent that further work is needed to understand why CMs feature 
in some cases serious biases on the integrated global budgets and disagree 
quantitatively quite substantially - while agreement on qualitative proper- 
ties is rather satisfactory - on the representation of a fundamental proper- 
ties such as the meridional enthalpy transport. Whereas in the atmosphere 
the agreement in the meridional atmospheric transport is reasonable, in the 
ocean rather different pictures emerge from the various CMs, especially in the 
southern hemisphere. In terms of response to CO2 increases, whereas we have 
relatively good agreement in the CMs response for the atmospheric transport, 
it is apparent that the ocean response requires a detailed re-examination. 
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The consideration of multi-model ensemble means has often been used as 
a way to partially circumvent the problem of dealing with CMs' discrepancies 
(see, e.g., [17, 59]). The rationale for this choice is that it is often found that 
such means provide enhanced skill with respect to randomly chosen specific 
CMs [13], even if this is far from being always the case [36]. Even if various 
statistical techniques have been proposed in place of the the usual arithmetic 
average [39, 40, 56] with encouraging empirical results, it must be emphasized 
that there is no rigorous statistical ground for combining data from different 
CMs [28] (outputs of different CMs are not samples of the same probability 
space, no large numbers law applies [30]), so that a detailed understanding 
of CMs' uncertainties and discrepancies is of the greatest urgency. 

In general, advances in numerical schemes as well as in the representa- 
tion of physical processes are needed for providing a well-defined and self- 
consistent representation of climate as a non-equilibrium thermodynamical 
system when long time scales and large spatial scales are considered. Fur- 
thermore, additional work is then needed in the definition of suitable metrics 
for the investigation of the global thermodynamical properties of the CS 
and in testing the performances of state-of-the-art and future generations 
of CMs. Along these lines, and stimulated by recent extensive scholarly re- 
views [43, 68], the authors are preparing a companion publication dealing 
with the analysis of the entropy production [43, 29] of PCMDI/CMIP3 CMs. 
Recently, [32] have shown that it is possible to estimate the total entropy 
production of the CS just by analyzing the TOA and surface longwave and 
shortwave radiative fluxes and to separate the contribution due to vertical, 
small scale processes from those due to horizontal, large scale processes. 

We suggest that the PCMDI/CMIP5 initiative, which will specifically 
take care of providing the tools necessary for auditing the present and future 
generation of climate and Earth system models, should take in greater care 
the provision of data referring to simulations by CMs under steady state 
with various choice of basic parameters, such as the solar constant (e.g., ± 
2% of the present one), the C0 2 (e.g., 275 ppm, 550 ppm, 1000 ppm), and 
should make it possible for research groups the access to data necessary for 
the analysis of the energy and entropy budgets of the planet. While these 
simulation and datasets might seem unnecessary and not of direct impact 
for the assessment of the actual and projected climate change, they would 
allow for a much more detailed understanding of the performance of CMs. 
Of course, it is crucial for matters of self-consistency that the energy bias 
of a CM is as small as possible in all steady state scenarios. It is apparent 
that climate-dependent energy biases may contribute to increasing the un- 
certainty in climate sensitivity, and we may conjecture that this may play a 
role in the relatively large discrepancies among state-of-the-art CMs on the 
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definition of this crucial quantity (still spanning between 2 and 4.5 C°), and, 
correspondingly, of the change in the globally averaged surface temperature 
resulting from historical and projected CO2 concentration increase [51]. 

We conclude by noting that most of the results presented here could also 
be applied for studying the climates, via models and observational data, of 
celestial bodies such as planets and satellites. This is a rather promising and 
exciting perspective, because, apart from the everlasting interest in the solar 
planets and satellites, we live in the dawn of the era of extra-solar planets 
and satellites science. Therefore, what reviewed here as well as the ongoing 
studies on the entropy production and efficiency of the Earth system may be 
of help for studying the thermodynamics of the atmosphere of celestial bod- 
ies. Thus, it is encouraging to observe that various models belonging to the 
PLASIM family have already been adapted to study the atmospheres of Titan 
[15] and Mars [52], and that intercomparison projects on the modeling of the 
Venusian atmosphere (see http : //www. issbern. ch/workshops/venusclimate 
and Martian atmosphere (see http://www.atm.ox.ac.uk/user/ 
newmanc/workshop/intercomp.html) are ongoing. 

The authors wish to thank M. Ambaum, E. Cartman, D. Battisti, K. 
Fraedrich, D. Frierson, J. Gregory, F. Lunkeit, S. Pascale, A. Speranza, and 
R. Tailleux for useful discussions over these topics. 
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Figure 1: First Panel: Idealized time dependence of the TO A global energy 
budget as given by a CM which has a self-consistent representation of energy 
exchange processes (solid line) versus an energy biased CM. Black, red, and 
blue lines refer to the initially unperturbed runs, the forced runs, and the 
runs performed with newly established fixed parameters, respectively. Second 
Panel: Smoothed (5-year filter) yearly time series of the GFDL2.1 model (CM 
9) TOA energy budget. The black line refers to the pre-industrial run. The 
red line refers to the XX century simulation (started from year 1 of the pre- 
industrial run) and to the XXI century portion of the SRESA1B simulations 
(started from the end of the XX century simulation). The blue line refers to 
the XXII and XXIII century SRESA1B simulation. 
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Figure 2: First Panel: Mean and standard deviation of globally averaged 
climate energy balance for all the models over 100 years long time periods, in 
the preindustrial scenario. Second Panel: Difference between the SRESA1B 
and the pre-industrial scenario. 
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Figure 3: First Panel: Mean and standard deviation of globally averaged 
atmosphere energy balance for all the models over 100 years long time pe- 
riods, in the preindustrial scenario. Second Panel: Difference between the 
SRESA1B and the pre-industrial scenario. The label a refers to data relative 
to the XXIII century in the SRESA1B scenario. 
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Figure 4: First Panel: Mean and standard deviation of globally averaged 
ocean energy balance for all the models over 100 years long time periods, in 
the preindustrial scenario. Second Panel: Difference between the SRESA1B 
and the pre-industrial scenario. The label a refers to data relative to the 
XXIII century in the SRESA1B scenario. 
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Figure 5: First Panel: Mean and standard deviation of globally averaged 
land energy balance for all the models over 100 years long time periods, in 
the preindustrial scenario. Second Panel: Difference between the SRESA1B 
and the pre-industrial scenario. The label a refers to data relative to the 
XXIII century in the SRESA1B scenario. 
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Figure 6: Mean total meridional enthalpy transport by the climate system 
(black lines), the atmosphere (red lines), and the ocean (blue lines) in the 
preindustrial scenario (first panel), and in SRESA1B scenario for the XXII 
century (second panel) for all models where suitable data are available. 
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Figure 7: Value and position of the peak of the total poleward meridional 
enthalpy transport in the preindustrial scenario for the northern hemisphere 
(first panel) and for the southern hemisphere (second panel). 
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Figure 8: Value and position of the peak of the poleward meridional en- 
thalpy atmospheric transport in the preindustrial scenario for the northern 
hemisphere (first panel) and for the southern hemisphere (second panel). 
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Figure 9: Value and position of the peak of the poleward meridional enthalpy 
ocean transport in the preindustrial scenario for the northern hemisphere 
(first panel) and for the southern hemisphere (second panel). 
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Figure 10: Difference in the peak value of the poleward meridional enthalpy 
transport for both hemispheres between the SRESA1B and the pre-industrial 
simulations. Data are provided for the total transport (first panel), the 
atmospheric transport (second panel), and the ocean transport (third panel). 
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Figure 11: Difference in the value and position of the peak of the atmospheric 
poleward meridional enthalpy transport between the SRESA1B and the pre- 
industrial simulations for the northern hemisphere (first panel) and southern 
hemisphere (second panel). 
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Figure 12: Difference in the value and position of the peak of the ocean 
poleward meridional enthalpy transport between the SRESA1B and the pre- 
industrial simulations for the northern hemisphere (first panel) and for the 
southern hemisphere (second panel). 



55 



Table 1: Models considered in this work, institutions and data availability in each considered time period. Each 
model is labeled with a number. 





Model 


Institution 


PI 


SRES A1B 
2101-2200 


SRES A1B 
2201-2300 


1 


BCCR-BCM2.0 


Bjerknes Center, Norway 


yes 


no data 


no data 


2 


CGCM3.1 (T47) 


CCCma, Canada 


flux adjustment 


flux adjustment 


flux adjustment 


3 


CGCM3.1 (T63) 




flux adjustment 


flux adjustment 


flux adjustment 


4 


CNRM-CM3 


Met?o France, France 


yes 


yes 


2101-2299 


5 


CSIRO-Mk3.0 


CSIRO, Australia 


80 years, no surf 


no data 


no data 


6 


CSIRO-Mk3.5 




yes 


no data 


no data 


7 


FGOALS-gl.O 


LASG, China 


yes 


yes 


no data 


8 


GFDL-CM2.0 


GFDL, USA 


yes 


yes 


yes 


9 


GFDL-CM2.1 




yes 


yes 


yes 


10 


GISS-AOM 


NASA-GISS, USA 


yes 


no data 


no data 


11 


GISS-EH 




no surf 


no data 


no data 


12 


GISS-ER 




no surf 


no surf 


2201-2297, no surf 


13 


HADCM3 


TT 11 /"I j T T t y 

Hadiey Center, UK 


yes 


2108-2199 


no data 


14 


HADGEM 




99 years 


2101-2199 


no data 


15 


INM-CM3.0 


T j C TVT ~K JT j 1 7~1 

Inst, oi Mum. Math, Russia 


yes 


yes 


i j 

no data 


16 


IPSL-CM4 


IPSL, France 


yes 


yes 


2201-2231 


17 


MIROC3.2 (hires) 


CCSR/NIES/FRCGC, Japan 


yes 


no data 


no data 


18 


MIROC3.2 (medres) 




yes 


yes 


yes 


19 


ECHO-G 


MIUB/METRI/M&D, Germany /Korea 


flux adjustment 


flux adjustment 


flux adjustment 


20 


ECHAM5/MPI-OM 


Max Planck Inst., Germany 


yes 


no surf 


no surf 


21 


MRI-CGCM2 


Meteorological Research Inst., Japan 


flux adjustment 


flux adjustment 


flux adjustment 


22 


NCAR CCSM 


NCAR, USA 


no surf 


no surf 


no surf 


23 


NCAR PCM 




no data 


2101-2199, no surf 


no data 



Table 2: Correlation between annual maxima of atmospheric and oceanic 
poleward transport in each hemisphere. The 95% confidence interval width 
is ±0.20 in all cases. PI and SRESA1B (2101-2200) scenarios. 
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4 


-0.32 


-0.29 


-0.46 


-0.42 


c 


6 


-0.36 


-0.26 






7 


-0.48 


-0.20 
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